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Abstract 



The open critical XXZ spin chain with a general right boundary and a trivial 
diagonal left boundary is considered. Within this framework we propose a simple 
computation of the exact generic boundary S'-matrix (with diagonal and non-diagonal 
entries), starting from the 'bare' Bethe ansatz equations. Our results as anticipated 
coincide with the ones obtained by Ghoshal and Zamolodchikov, after assuming suitable 
identifications of the bulk and boundary parameters. 
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1 Introduction 



The derivation of exact bulk and boundary S-matrices in the context of 2-dimensional in- 
tegrable models is a fundamental problem that has attracted considerable attention during 
the last decades (see e.g. p] [6]). The main aim of the present investigation is the compu- 
tation of the exact generic -with diagonal and non-diagonal entries- boundary S-matrix in 
the framework of the open XXZ spin chain. More precisely, we consider the open XXZ spin 
chain with the left boundary to be trivial (K + oc I), while the right boundary is associated 
to a full iT-matrix J5[ [7] . The corresponding Hamiltonian reads as: 

n = + Oi<>i+i + cosh(z/i) ofo^j) - j cosh(2/i) - Smh ^ <T^ 

1=1 

sinh(m) cosh(in£) Ksinh(m) / . . \ 

where cr a: ' 2/ ' 2: are the 2x2 Pauli matrices, and the boundary parameters £, k, 8, are the free 
parameters of the generic K-mathx and are associated to some magnetic field applied at the 
boundaries of the chain. Notice that we focus here on the critical regime |e JM | = 1. 

In general quantum spin chains are one dimensional statistical systems displaying particle- 
like excitations 'holes' which may scatter among each other [2] or with the boundary [SJ, IS] 
in the presence of non trivial boundary magnetic fields. The boundary S-matrix within this 
framework describes exactly the reflection of the the particle-like excitation with the bound- 
ary magnetic field. Note that computations concerning both bulk and boundary S-matrices 
are always valid in the thermodynamic limit of the spin chain iV — > oo [2J m [TTJ [SJ M, ED] • 
In integrable field theories on the other hand (see e.g. [5]) the physical boundary S'-matrix 
describes the reflection of a solitonic excitation with the boundary, when the system is con- 
sidered on the half line [5] (see also [121 IS] for relevant studies on the sine-Gordon model 
on the finite interval). It is crucial to note that the critical XXZ model may be thought 
of as the discrete analogue of the sine-Gordon model, see for instance [HI [15[ [16] for a de- 
tailed derivation of this correspondence. It was also shown that the bulk [171 [10] and the 
diagonal boundary S'-matrices [HI [10l [9] of the critical XXZ chain coincide with those of 
the sine-Gordon model after considering suitable identifications of both bulk and bound- 
ary parameters. It is therefore naturally anticipated that the generic boundary S-matrices 
should also coincide. Such a coincidence is also natural from the algebraic point of view since 
both models are ruled by the same algebra defined by the reflection equation. Recall also 
that boundary S-matrices are solutions of the reflection equation up to an overall physical 
factor, which in the spin chain frame may be exactly computed by means of Bethe ansatz 
techniques. 
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There exist several schemes to derive physical boundary S'-matrices, such as the bootstrap 
method developed in [HE], the 'non-linear integral equation' (NLIE) technique (see e.g. [EJ 
[12], [20]), and the 'physical' Bethe ansatz formulation [TTJ, [TH] - However, the most direct means 
is arguably provided by the 'bare' Bethe ansatz approach (see for instance [H [HI El [TO]). 
In this study the derivation of the generic boundary iS-matrix is based on the 'bare' Bethe 
equations [2"Tj |2"2"] . and this is the first time that a direct computation starting from the 
microscopic Bethe ansatz equations is achieved for the generic reflection matrix. 

The crucial observation is that in the special case under consideration the Bethe equations 
reduce to an elegant and quite familiar form after assuming suitable boundary parametriza- 
tions parallel to the ones of [5]. Specifically, the Bethe ansatz equations acquire a form 
similar to the one of the open XXZ spin chain with two diagonal boundaries (see also [23] 
for relevant comments). We are hence able to compute the reflection amplitudes in a simple 
fashion analogous to the diagonal case, and this is actually one of the main advantages of 
the approach adopted here. It is worth pointing out that in the case where two generic 
boundaries are implemented the entailed Bethe ansatz equations are rather involved [2U [25] 
without really offering any additional information as far as the boundary scattering is con- 
cerned. Note that in [25] the spectrum in the most generic case with general free bulk and 
boundary parameters is obtained. Given the complexity of the Bethe ansatz type equations 
in the most general case it is clear that our choice of simpler boundary conditions is quite 
practical rendering the relevant computations considerably more tractable. Moreover, in 
the particular case we assume here there exists a non-local boundary conserved quantity, 
at first order, contrary to the generic case. Comments regarding its role on the spectrum 
and Bethe ansatz equations, and on its relevance to S z are also presented. Our results as 
expected coincide with the ones obtained by Ghoshal-Zamolodchikov for the sine-Gordon 
model with 'free' boundary conditions [5j. Relevant results from the NLIE point of view 
were also derived in [20] . 



2 The open XXZ spin chain 



We recall the i?-matrix associated to the the spin 
Yang-Baxter equation [26]. 



| XXZ model, which is a solution of the 



Rn(^i — A 2 ) i?i3(Ai) #23(^2) = -^23(^2) -^13(^1) Rui^i — A 2 ). 



(2.1) 



The i?-matrix has the following form: 




(2.2) 
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it is convenient for our purposes here to set // = ~. 

The transfer matrix of the open spin chain is defined as [2 

t(X) = tr {M K + (X) T(A) K~(X) f (A)} 

T(X) = R 0N (X) ...R 01 (X), f(X) = R w (X)...R N0 (X) (2.3) 

where M = diag(q, g _1 ), q = e lfl and K ± are solutions of the reflection equation |28j : 

i? 12 (Ax - A 2 ) K 1 {X 1 ) R 21 (X 1 + A 2 ) K 2 (X 2 ) = K 2 (X 2 ) R 12 (X 1 + A 2 ) K 1 (X 1 ) R 21 (X ± - A 2 )(2.4) 

The general solution K(X) is a 2 x 2 matrix with entries [5] [7]: 

ifn(A) = sinh[/i(-A + zO]e M , K 22 {X) = sinh[^(A + *0]e~ MA 

K 12 (X) = Kg sinh(2//A), K 2 i(A) = sinh(2//A). (2.5) 

We shall henceforth consider = I, and K~(X) = K(X) defined in (12.51) . This particular 
choice of boundary conditions is compatible with certain constraints imposed upon the left 
and right boundary parameters [22], [21], [29], [30] . We are eventually left with three free bound- 
ary parameters associated only to the right boundary. As a matter of fact the parameter 9 
may be removed via a simple gauge transformation, as also happens in [S], where a X-matrix 
with 8 = is assumed for simplicity, but without loss of generality. Using the fact that the 
quantity (T K~ T) satisfies the reflection equation one can show that the transfer matrix 
(I2.3P provides a family of commuting operators [27j : 

t(A), t(A')j = 0. (2.6) 

Evaluating the eigenvalues of the i^-matrix will be particularly useful for both adopt- 
ing an appropriate parametrization for our Bethe ansatz equations as well as for comparing 
effectively with the Ghoshal-Zamolodchikov result. Before writing down the .fT-matrix eigen- 
values it will be useful to introduce some notation: 

— — = *coshM/r + 7~)], — = *cosh(z/<), p ± = _(/T+ 7 -±C) (2.7) 
2k 2k 2 

compare also with the parametrization used e.g. in [20] [21] [22] [29] [30]. Such a parametriza- 
tion is also quite natural from the point of view of boundary Temperley-Lieb algebras [31] . 
The latter formulas (12. 7p lead to the following relations among the boundary parameters: 

CO Sll ( 2 Lie ) 1 

= cosh(ifip + ) cosh(z^j9~), cosh 2 (z//p + ) + cosh 2 (z/ip~) = 1 (2.8) 

2%k Ak 2 

(c.f. see similar parametrization in [5], and section 3). 'Renormalized' physical boundary 
parameters will be ultimately identified with the boundary parameters of |5], as will become 
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transparent in section 3. The parameter 9 is hidden is the sum (3~ + 7" (see e.g. [22| 29J), 
but in any case we shall consider it henceforth to be zero for simplicity. After diagonalizing 
the i^-matrix we obtain the two eigenvalues: 

£l 2 (A) = -2in sinh[/i(A ± ip + )] sinh[/x(A ± ip~)\. (2.9) 

It is worth stressing that the Bethe ansatz equations provide essentially the eigenvalues of 
the physical boundary S-matrix -with 'renormalized' boundary parameters. Nevertheless, 
the structure of the boundary S'-matrix is a priori known due to the fact that is a generic 
solution of the reflection equation (12. 4p . The objective now is to obtain the overall physical 
factor in front of the boundary S'-matrix; this will be achieved in the subsequent sections by 
means of the Bethe ansatz approach. 



2.1 Bethe ansatz equations 

The spectrum and Bethe ansatz equations for the spin | XXZ chain with general boundary 
conditions were derived in [221 EE] , whereas in [291 EH E2] the spectrum and Bethe ansatz 
equations obtained for various representations. We shall focus here on the spin | case, the 
spectrum in this case is given by: 

A(A) = A-(0|A)A 7 ( P ±|A) sinh- M A + i)] f\ + ^ f^u Til 

w 1 K j=i simi [M^ + \ + l )] smh[/i(A - Xj)] 

+ Ktmxw *v w u^B^W^W 1 (210) 

where Kf±, for the particular choice of boundary conditions are given by (see also Appendix 
in [29] for the explicit definitions of K^ 4 ) 



Kiip^X) = -2ik e MA sinh[y[x(A - ip~)] sinhf^A - ip + )], 
K 4 (p ± \X) = -2in e MA sinh[//(A + ip~ + i)] sinh[/i(A + ip + + i)] 



sinh(i//) 



Jtf (0|A) = e-* Sin !;i 2 ^ + ^ , Kt(0\X) = e~>* „ . (2.11) 

1V 1 ' sinh[/x(2A + i)} ' 4V 1 ' sinh[//(2A + i)} K ' 

The Bethe ansatz equations arise as necessary constraints such that certain 'unwanted' 
terms appearing in the eigenvalue expression are vanishing. They guarantee also analyticity 
of the spectrum, and are written in the familiar form: 

M 

e 2p-+i( A 'i/ i ) e 2p++i( A ^) 9i{h\IJ) el N+1 (Xi;n) = - JJ e 2 (Xi - Xj; fi) e 2 (A i + A i ; y u) , (2.12) 
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where we define 



in x 



sinh[MA + f )] cosh[/i(A+ 
e " (A ^ ) = sinh[MA-f)] ' ^ (A ^ ) = coshKA-f)] - (2 ' 13) 

The similarity of the Bethe ansatz equations (I2.12p with the Bethe equations of the case 
with two non-trivial diagonal boundaries is indeed noticeable (see e.g. [HI El EZJ)- This is a 
crucial observation, which will considerably simplify the derivation of the generic boundary 
^-matrix. 

In [221 [291 EO] the spectrum and Bethe ansatz equations were derived starting from a 
particular reference state (the analogue of 'spin up' state), however it was shown in |33j that 
a second reference state exists ('spin down'), providing another set of Bethe ansatz equations 
similar to (I2.12p . but with — > —p ± , guaranteeing also the completeness of the spectrum. 
For a detailed analysis on the completeness of the spectrum using both reference states we 
refer the interested reader e.g. to [2U [331 E2]- It should be stressed that the existence of 
the two sets of Bethe ansatz equations together with the similarity of f l 2 . 1 2 f) with the Bethe 
equations in the purely diagonal case are the key elements in deriving the exact boundary 
.S-matrix. 

Interestingly in the case we are considering here, where the left boundary is trivial, there 
exist a conserved quantity, which is somehow the analogue of the spin S z of the diagonal case 
|27j . The integer M appearing in the Bethe ansatz equations is actually associated to the 
spectrum of the conserved quantity (see e.g. [29J). Let us briefly recall the structure of the 
boundary conserved quantity (boundary non-local charge), which is defined as [31 



Q(N) = q -l 2+ e K (N) E (N) + q l-e K (N) F{ N) e _^( K (N) ) 2_ (2 M) 

2k sinh(z/i) 

K( N \ EW and are N coproducts of the quantum algebra U q (sl2) [37] i.e. 

N N 

K {N) = (g) K n , X {N) = K- 1 <g> . . . R- 1 <g> X ®K...®K, 



n=l n=l 



n position 



X e {E, F}. (2.15) 
For the spin | representation in particular [K — > q 2 ^ , E — > cr + , F — > a~) we have: 

Q (N) = ~ ■ uV \ cosh^(/5- + 7 " - 2S)}. (2.16) 
smh(zyu) 

The one-site Q^ N ' operator becomes a 2 x 2 matrix and the eigenvalues of S are ±~. For 
N = 2, S has 0, ±1 as eigenvalues, with being a doubly degenerate eigenvalue, and 
so on. It is thus quite clear that the operator S behaves similarly to S z . Details on the 
diagonalization of the operator Q( N > for the spin | representation can be found in [36J. 
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From the asymptotic behaviour of the open transfer matrix (12. 3p and the asymptotics 
of the spectrum (12.101) the following fundamental formula emerges (for more details on this 
matter we refer the reader to |29j): 

Q( N ) = - . J cosh^/T + 7- + 2M - N)}. (2.17) 

the subscript e in the latter expression stands for the eigenvalue. In this case M is not fixed, 
as opposed to the generic case , but is associated to the spectrum of the operator Q^ N \ From 
equation (I2.17P the upper (lower) bounds of M (integer) are identified from the spectrum of 
the non-local operator Q^ N \ indeed: 

N 

M=--S £ . (2.18) 

The similarity of the latter relation with the one appearing in the case of diagonal boundaries 
is noticeable; indeed in the diagonal case S is simply replaced by S z in (12.181) . In fact, the 
symmetry in this case is effectively U(l) resembling the case of purely diagonal boundaries 
(see also [38]). In general the identification of the spectrum of the boundary non-local charge 
Q( N ^ is an intriguing problem and particular cases have been analyzed in [361 138] . 

Finally, comparing the eigenvalues of (12.161) for an one particle state (N odd) -the eigen- 
values for S are then ±|- with (12.171) we conclude that: 

N - 1 

M=^-. (2.19) 

This corresponds to the state with S eigenvalue |, the other state with eigenvalue — | -above 
the equator- may be obtained via a 'duality' transformation on the boundary parameters, 
as it happens in the diagonal case [27J; this will be demonstrated however in the subsequent 
section. 



3 The boundary /S-matrix 

As already mentioned in the introduction the physical boundary S matrix describes the 
reflection of a particle-like excitation (kink), displayed by the XXZ spin chain, with the 
boundary. Our main objective here is the computation of the overall physical factor in front 
of the boundary S'-matrix, which contains significant information regarding the existence of 
boundary bound states. 

The physical boundary S'-matrices are denoted henceforth K . We define the boundary 
^-matrices by the quantization condition [TT], [S] 
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A is the rapidity of the 'hole' -particle-like excitation, and p(X) is the momentum of the hole, 
which will be formally defined shortly. Recall that both bulk and boundary scattering may 
be validly evaluated in the thermodynamic limit (N — > oo) [21 (8J [9j HQ] . In this limit as is 
well known the so called string hypothesis is valid [2]. Based on this hypothesis (see e.g. 
[2]) it was shown that the ground state for the antiferromagnetic spin chain, which we study 
here, consists of a 'sea' of real (1-string) Bethe roots, and a particle like excitation is simply 
a 'hole' in this uniformly distributed sea of real roots [21 HI E]. For the Bethe ansatz state 
with one hole of rapidity A in the sea, the counting function -obtained after taking the log 
of Bethe ansatz equations- is 

h{X) = ^\(2N+l)q 1 (X;fi)+r l (X;fi)-q 2p++1 (X;fi)-q 2p _+i(X;fi) 

M ^ 

fe(A - A, ; fi) + q 2 (X + X, ; fi)] j, (3.2) 



where 



q n {X;fi) = 7r + ilog[e n (A;/i)], r„(A; fi) = i log[#„(A; fi)]. (3.3) 



Note the striking similarity of the latter formula (13.21) with the one obtained in the case of 
two diagonal boundaries [HI E]- This is a key point, as already mentioned, allowing us to 
proceed with a simplified derivation of the boundary S-matrix. The only difference with 
the fully diagonal case is that now both terms depending on p are assigned to the right 
boundary, otherwise we proceed exactly as in the diagonal case (see e.g. [HIE])- 

Since we are considering the thermodynamic limit it is clear that all the sums in the ex- 
pressions above turn into integrals, and also we need a density to describe the corresponding 
states, see e.g. [21 El El ED] for details. The Fourier transform of the density er s (A) = jj^jp- 
describing the one-hole state is given by 

a s (uj) = 2e(uj) + — — — r(e +e ) 

N 1 + a 2 {uj ; fi) 



1 1 

+ 



1 + a 2 (uj ; fi) 



ai(u ; fi) + a 2 (u ; fi) + bt(uj ; fi) - a 2p - +1 (uj ; fi) - a 2p++1 (u ; fi) , 

(3.4) 



where we define the following Fourier transforms (recall that |e* M | = 1 and fi = ^): 

a n (u; fi) = s r ih{ ^ 21 0<n<2u (3.5) 

f , . sinh(^) / 

b n (u-fi) = - ginh ^^ 0<n<u (3.6) 
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and 



e w 



(3.7) 



l + a 2 (w;/i) 2 cosh (f)' 

e(A) corresponds also to the energy of the particle-like excitation. 

The boundary matrix K~, of the generic form (I2.5p . has two eigenvalues k 12 . The left 
boundary matrix is trivial 



K+(A) = fco(A)I. 



(3i 



From the density (13.41) and the quantization condition (13. ip we explicitly derive the quantities 
k , ki >2 . Indeed taking into account that 



6(A) 



1 dp{\) 



(3.9) 



2tt dX 

and comparing ( 13.4f) with the quantization condition ( 13. II) we obtain the first eigenvalue ki 
(see e.g. [H [9] for more details). It is worth noting that the momentum is not a conserved 
quantity anymore, however the momentum of the particle-like excitation may be formally 
defined via (J33D (see also [8, 9J). 

To obtain the second eigenvalue k 2 we apply the argument of [271 [9] for the diagonal case 
i.e. implement a 'duality' transformation on the boundary parameters such that p* 1 — > — p ± . 
This is actually the equivalent of deriving the Bethe ansatz equations starting from the 
second reference state ('spin down') [33l 132] . The explicit expression for the eigenvalue ki is 
given by: 



ki(A, p + , p ) = 
where we define 



-%K 



cosh 



7T^ 



V 



-IV 2^ 



2p H 



cosh 



7T 



V 



A 



fc(A, p + , p-) = k (X) h(X, P + ) h(X, p~) 



-{v-2p )) k(X,p + p ), 



(3.10) 



and 

fci(A, x) 



exp <^2 ^ — sinh ^2iu;A j 



sinh(^) sinh ((z/- 2)f) 



7r(— 2zk) 2 



cosh 



A-i(*-2s) 



sinh (2w) sinh ((u - 1)|) 

°° da> • , r\ sinh ((i/ — 2x — l)u 
- sinh ( IiujX ) — 

'o w 



exp | — 2 / — sinh ^2icjA j 



2 sinh ((i/ — l)u;) cosh(iu; 

(3.11 
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The latter expressions may be written as infinite products of T-functions, so the comparison 
with the results in [5] is easier 



fco(A) 



x 



r (^(-2zA + An + 3) + l) T (^(-2iA + An + 1)) 



n=0 I r (^(2zA + 4n + 3) + lj T (^(2zA + 4n + T 

T (^(2zA + 4n) + l) T (^(2zA + 4n + 4) 
T (^(-2zA + An) + l)v (^(-2^ + An + 4) 



(3.12) 



fci(A, x) 




x 



T (^(-iA + 2n - \{y - 2x)) + f) T C^(-iX + In + - 2ar)) + \ 



H I T (^(iA + 2n + 2 - \{y - 2x)) + |J T ( ^(iA + 2n + 2 + - 2x)) + \ 
r (^(iX + 2n + 1 - - 2a;)) + |) T (^(iX + 2n + 1 + \{v - 2x)) + \ 



(3.13) 



r s=i (~«A + 2n + 1 - f (i/ - 2x)) + | T ^(-«'A + 2n + 1 + §(v - 2x)) + | 



After implementing the 'duality' transformation on the boundary parameters p ± — > —p in 



Bethe ansatz equations ( 12.121) we obtain the second eigenvalue k 2 of the reflection matrix 
with: 



cosh 



k 2 (A, p + , p ] 

ki(A, p+, p-) cos h 



E T (X+i(i/-2p+))] cosh[^(A + |(i/-2p 



j^j A- f(i/-2pH 



cosh 



-MA 

I/— 1 ' 



- 2p- 



(3.14) 



compare also the latter relation with (12.91) ; the appearance of renormalized parameters be- 
comes now apparent. Notice that the p dependent term, expressed as a product of two k\ 
functions (13.131) . is 'double' compared to the diagonal case. In the diagonal limit only one 
of the functions &i(A, p^) survives, while the other one becomes unit (see also relevant dis- 
cussion in [5]). The diagonal case corresponds to what is called 'fixed' boundary conditions 
in 0. 

Comparison between our findings and the results of [5] for 'free' boundary conditions, 
gives rise to the following identifications of bulk and boundary parameters. The Ghoshal- 
Zamolodchikov bulk coupling constant^ A is related to our coupling constant v by A = -^-r] 
the spectral parameter u in [5] is related to our variable A by u = —mX. Note that the 
original parameter fi = ^ appearing in the i?-matrix (I2.2p is now renormalized to jl = 



2 The coupling constant A is denned in [5] as A 
constant. 



1, where (3 is the familiar sine- Gordon bulk coupling 
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Recall also that the boundary parameters in [5] (rj, •&) and (k, £') (we denote £' the Ghoshal- 
Zamolodchikov ^-parameter to distinguish it from the bare parameter of (12.51) ) satisfy [5]: 

cos(r?) cosh(i9) = — — cos(£'), cos 2 (??) + cosh 2 ^) = 1 + — , (3.15) 
k k 2 

which are similar to the constraints (I2.8P among the 'bare' boundary parameters. Finally, 
the following identifications among the boundary parameters are valid: 

m{v-2p+) rc(u-2p-) 7TQ/-2Q 

2(i/-l) ' 71 = 2(u-l) ' C= ^l7' * = ~ 2 *« ( 3 - 16 ) 

see also relevant formulas in [20]. With this we conclude our derivation of the generic physical 
boundary S-matrix, which naturally coincides with the general reflection matrix found in [5] 
associated to 'free' boundary conditions. 



4 Discussion 

Let us summarize the main findings of this investigation. Our main objective was the 
derivation of the exact generic boundary S'-matrix for the open XXZ chain. This was achieved 
by means of the 'bare' Bethe ansatz approach. More precisely, we considered a particular 
case of the open spin chain, with a generic right boundary and a trivial left one. Then 
assuming appropriate boundary parametrizations we were able to write down a simple and 
familiar form of the Bethe ansatz equations, similar to the Bethe equations of the XXZ chain 
with two diagonal boundaries. 

The simple form of the Bethe ansatz equations has been a crucial point in our analysis. 
This together with the existence of a second set of Bethe ansatz equations, obtained via 
a duality transformation on the boundary parameters, facilitated the derivation of the two 
eigenvalues of the generic physical boundary S'-matrix. Moreover, we were able to associate 
the spectrum of the boundary conserved quantity Q^ N ' with the spectrum of the transfer 
matrix, and we showed that plays essentially a role analogous to S z in the fully diagonal 
case. We identified the overall physical factor in front of the boundary S-matrix of the form 
(12. 5p , but with renormalized boundary parameters (13.161) , which are basically the parameters 
used in |5J. Thus the reflection matrix derived in [5] was fully recovered. 

Similar results may be deduced in the attractive regime where breathers (bound states) 
are also present [10J. In this case the computation of the generic breather boundary S-matrix 
goes along the same lines as in |10| . The n breather boundary S-matrix is then given by: 

Sl n \\, p+, p-) = (A) s{ n) (A, p + ) s{ n) (A, p-) (4.1) 
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where explicit formulas for Sq (A), Si (A, x) are presented in [10]. Expression (14. ip is 
analogous to the n th breather reflection matrix in the sine-Gordon model obtained in [39J for 
'free' boundary conditions. In the diagonal limit one of the terms in ( 14.1ft becomes unit, 
exactly as in the solitonic reflection matrix, and the result of [10J is recovered corresponding 
to the 'fixed' boundary conditions in |39j. It is clear that such a computation immediately 
provides the reflection matrix for the fundamental particle in sinh-Gordon. We recall that 
the lightest breather (n = 1) in sine-Gordon corresponds essentially to the fundamental 
particle of the sinh-Gordon model, provided that the coupling constant (3 — > i/3. 

Another problem which may be treated in the same spirit is the derivation of the generic 
reflection matrix for the XXZ chain in the non-critical regime. In this case the boundary 
S'-matrix will be expressed in terms of I^-functions, see e.g [9j HQ] for diagonal boundaries 
only. Finally an interesting direction to pursue is the computation of generic boundary S- 
matrices in the context of higher spin open XXZ chains. Specifically, the spin-1 case [T71 14"T] 
is of particular significance given its relation to the super-symmetric sine-Gordon model [42J. 
We hope to address these issues in forthcoming publications. 

Acknowledgments: This work was supported by INFN, Bologna section, through grant 
T012. 
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